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I. INTRODUCTION 



Recently we reported [1] on the existence of different intermediate Hamiltonians in type 
A A/"-fold supersymmetry (SUSY) [2]. Apart from the connectivity to a wide range of 
solvable and quasi-solvable potentials including not only the well-known sl(2) Lie-algebraic 
models [3] but also the new extended class of completely solvable rational ones [4], type 
A A/"-fold SUSY is of independent interest because of the rich mathematical structure it 
provides (for a review, see Ref. [5]). 

The main motivation of our study came from the issue of non-uniqueness of factorizing 
operators as a consequence of the underlying GL(2, C) symmetry. Type A Af-iold super- 
charge admits of a one-parameter family of factorizations that is expressible as a product 
of M first-order linear differential operators due to the aforementioned symmetry [6] . This 
in turn implies that a type A A/"-fold SUSY can have different intermediate Hamiltonians 
corresponding to different factorizations. However the existence of intermediate Hamiltoni- 
ans is not guaranteed in general [2]. In Ref. [1], we investigated under what conditions type 
A TV-fold SUSY systems can admit intermediate Hamiltonians and that how many sets of 
such Hamiltonians are plausible for the specific case of Af = 2. We then concluded that the 
number of admissible intermediate Hamiltonians would be a more suitable index than the 
concept of reducibility introduced in Refs. [7, 8] to characterize higher-order intertwining 
operators. Furthermore, we found that it naturally leads to a realization of second-order 
parasupersymmetry (paraSUSY) [9] and generalized 2-fold superalgebra [10]. 

In this article, we pursue our investigations further and focus on the Af = 3 case. In 
comparison with the previous Af = 2 case, there appear mainly two novel features in the 
M = 3 case. The one is the fact that the functional types of type A A/"-fold SUSY potentials 
for Af > 3 are strictly restricted to at most elliptic functions due to an additional constraint. 
But due to this constraint, all the type A A/"-fold SUSY models for Af > 3 were completely 
classified in Ref. [6]. We find that we can also classify entirely all the type A 3-fold SUSY 
models which have intermediate Hamiltonians. 

The other novel feature of Af = 3 is concerned with the fact that there are two different 
intermediate positions in each factorized 3-fold supercharge. Due to the latter fact, one can 
consider different patterns in the existence of intermediate Hamiltonians. That is, in cer- 
tain cases systems admit intermediate Hamiltonians at both the two intermediate positions 
while in other cases systems admit them at only one of the two positions. Explorations 
of the underlying conditions for each case reveal that for the former classes a system with 
intermediate Hamiltonians at the two positions turns out to be always solvable and even 
shape invariant. On the other hand, for the latter classes a system with an intermediate 
Hamiltonian at only one position is led to quasi-solvability only. 

Thanks to the aforementioned different patterns in the Af = 3 case, we further find 
intriguing and quite rich structure which does not exist in the Af = 2 case when we consider 
the existence of more than one sets of intermediate Hamiltonians. For instance, we find 
that there are systems which have intermediate Hamiltonians at the two positions in one 
factorization but has only one at one of the positions in another factorization. Throughout 
the analyses of such systems, we realize that in contrast to the Af = 2 case we must consider 
not only the number of admissible intermediate Hamiltonians but also the variety in the 
existence of them to characterize A/th-order intertwining operators for the Af > 3 cases. 

We also investigate a parafermionic formulation of such systems and realizations of third- 
order paraSUSY [11, 12] and generalized 3-fold superalgebra [13]. We find not only that 
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such realizations are indeed possible but also that variations of the latter superalgebra hold 
according to the different patterns in the existence of intermediate Hamiltonians. 

The article is organized as follows. In Section II, we review type A 3-fold SUSY and 
discuss some of its salient features. Then, we introduce three classes according to three 
different patterns in the existence of intermediate Hamiltonians. In Section III, we work out 
explicitly the necessary and sufficient conditions for each of the three classes. In Section IV, 
we classify exhaustively different type A 3-fold SUSY potentials with one set of intermediate 
Hamiltonians. In particular, we find an intimate relation between the patterns and the degree 
of solvability. That is, all the systems which have intermediate Hamiltonians at the two 
positions consist of shape invariant potentials. In contrast, those which have one at only one 
position comprise sextic anharmonic oscillators, deformed Poschl- Teller or Scarf potentials, 
and one-body elliptic Inozemtsev potentials, all of which are quasi-solvable. In Section V, we 
discuss and study in full detail the necessary and sufficient conditions for a system to have 
different sets of intermediate Hamiltonians in various patterns. In Section VI, we then give 
the complete classification of the type A 3-fold SUSY models which admit simultaneously 
more than one sets of intermediate Hamiltonians. In Section VII, we show that a system with 
intermediate Hamiltonians at the two positions always admits a realization of third-order 
paraSUSY. In addition, each system belonging to one of the three classes always admits a 
realization of variant generalized 3-fold superalgebras. In the final section, we summarize 
and discuss various aspects of the obtained results. 

II. TYPE A 3-FOLD SUPERSYMMETRY 

A type A 3-fold SUSY system is given by 

H± = -\^i + \ W ^? ~ \ ( 2E '(x) ~ E ^f) ~ R± \ W '^ ( 2 - la ) 

^3 = -^31-^32-^33) -^3 + = — (-^3 ) = -^33-^32-^31' (2.1b) 

where the superscript T denotes the transposition in the rc-space and P^ (i = 1,2,3) are 
defined by 

p± = T d + W-E, P± = T d + W, P± = Td + W + E. (2.2) 
In the expanded form, the type A 3-fold supercharge component P 3 ~ reads 
p- = d 3 + 3Wd 2 + (3W + 3W 2 + 2E' - E 2 ) d 

+ W" + 3WW + W 3 + (2E' - E 2 ) W + - (2E' — E 2 )' . (2.3) 

The functions E(x) and W(x) are not arbitrary but are connected with a fourth-degree 
polynomial A(z) and a second-degree polynomial Q(z) 

A(z) = a^z 4 + a^z 3 + a 2 z 2 + a\Z + a , (2.4a) 
Q{z) = b 2 z 2 + b ± z + b , (2.4b) 

through the following relations 

z"(x) = E(x)z'(x), 2A(z) = z\x) 2 , Q{z) = -z'(x)W(x). (2.5) 
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From these relations, we obtain in particular 

W\x) + E(x)W(x) = -Q'(z) 



E'(x) + E(x) 



A"(z) 



(2.6) 



which will play important roles in later analyses. The restriction on the form of A(z) arises 
in type A A/"-fold SUSY for all M > 3 and is absent for Af = 2. This strongly limits the 
possible forms of potentials for Af > 3. 

Due to the underlying algebraic structure of type A A/"-fold SUSY systems, they are 
invariant under the linear projective transformations of GL{2, C) defined by 



aw + j3 



(a, 07, <SeC, A = a5 - fa ^0). 



(2.7) 



The polynomials A(z) and Q(z) are then transformed under the GL(2, C) transformations 



as 



A(z) ^ A(w) = A~ 2 ( 7 u> + 5) 4 A(z) 
Q{z) ^ Q{w) = A^^w + 5fQ{z) 



" -yw+5 k=Q 



_ qw-\-{3 



h k w k , 



E 

k=0 
2 



(2.8) 
(2.9) 



where the new coefficients hi (i = 0, . . . , 4) and bi (i = 0, 1, 2) are respectively given by 



/ 0, 4 \ 

a 3 

a x 
\h / 



( 



a 



a 3 7 



4a 3 /3 a 2 (a<5 + 3/37) 

6a 2 /3 2 3a/3(a<5 + /3 7 ) 

4a/3 3 P 2 {3a5 + fa) 
fa (3 3 5 



a 2 7 2 



cry 7 

2«7(a5 + /?7) 7 2 (3a5 + /?7) 47 3 <5 

a 2 5 2 + 4a/3 7 <5 + /3 2 7 2 3 7 <5(a<5 + /3 7 ) 6 7 2 <5 2 

2/3c>(a<5 + /3 7 ) 5 2 (a5 + 3/3 7 ) 4 7 5 3 
/3 2 5 2 /3<5 3 



4 \ 



54 I 





( CL4 \ 




a 3 




a 2 







(2.10) 



and 



cry 7 
a/3 a<5 + fa 2^5 
(3 2 (35 5 2 

From the formulas (2.5), the induced transformations of functions E(x) and W(x) read 

2^z'(x) 



(2.11) 



W(x) ^ lU(x) = W(a;), E(x) ^ £(x) = £(x) 



(2.12) 



^z[x) — a 

The latter transformation rule in particular implies the invariance of the following function: 

2E'{x) - E(x) 2 = 2E'(x) - E(x) 2 . (2.13) 
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The GL(2, C) invariance of any type A 3-fold SUSY system is now manifest since both the 
pair of Hamiltonians H ± in (2.1a) and the operator P 3 in (2.3) only depend on W and 
2E' — E 2 which are both invariant under any GL(2, C) transformations. 

The superHamiltonian H 3 and the type A 3-fold supercharges Q 3 introduced with the 
ordinary fermionic variables V ;± as 

H 3 = H-^-^ + + H + ^j + 4j-, Qf = P 3 ^ ± , (2.14) 

satisfy the type A 3-fold superalgebra [6] 

[Qf,H 3 ] ={Q-,Q+}=0, (2.15a) 

{Qs, Qt} = 8(H 3 + Rf~l WA] - 3D 2 [Q]) (H 3 + R) 

+ ^(HA]+9h,2[A,Q]), (2.15b) 

where i 2 [A], D 2 [Q], js[A], and Ii t2 [A,Q] are the absolute invariants composed of A(z) and 
Q(z) as the followings (cf., Ref. [6] and the references cited therein): 

i 2 [A] = 12a a 4 - 3aia 3 + a 2 2 , D 2 [Q] = 4b b 2 - ft 2 , 
2j 3 [A] = 12a^a 2 a^ — 27aoa 3 — 21a\a^ + 9aia 2 a 3 — 2a 2 , (2-16) 
h,2[A, Q] = 6a 4 5 2 - 3a 3 6 6i + 2a 2 b b 2 + a 2 b\ - Za^b x b 2 + 6a 6 2 2 . 

One of the most important aspects of the type A 3-fold SUSY system (2.1) is that its 
components satisfy the third-order intertwining relations: 

-^31-^32-^33-^ = -^31-^32-^33' ^33 ^32 ^31 ^ + = ^ ^33-^32-^31' (2-17) 

which are responsible for the first commutation relation in (2.15a). Before considering 
the central issues on the existence of intermediate Hamiltonians, we shall refer to another 
transformation property of the system. That is, the type A 3-fold SUSY system (2.1) is 
transformed under W(x) — > — W(x) as follows: 

H <H> H + , P 31 <H> —P 33 , P 32 -H- — P 32 , P 33 <->■ — P 3 \, P 3 <->■ — P 3 ■ (2-18) 

Hence, the system as a whole remains invariant under the transformation W(x) — > —W{x). 
In particular, it leaves the pair of intertwining relations (2.17) invariant. This transformation 
however does not belong to the GL(2, C) transformation since it changes the sign of W(x) 
while any GL(2, C) transformation does not change W(x) at all as shown in (2.12). 

For type A 3-fold SUSY systems, there are essentially three different patterns in the 
existence of intermediate Hamiltonians according to which we shall classify them as follows: 



Class (1,1): 



P 33 H- = H ll P 33 , P 3 \H il = H-P+, (2.19a) 
P 32 H il = H^P 32 , P+W l = H ll P 3 + 2) (2.19b) 
P 31 W l = H + P 31 , P+H+ = H^P 3 \. (2.19c) 
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Class (0,1): 



-^32-^33-^ — H^P 32 P 33 , P 33 P 32 H sl — H P 33 P 32 , (2.20a) 

P 31 W l = H + P 31 , P 3 \H + = H^P 3 \. (2.20b) 

Class (1,0): 

P 33 H- = H ll P 33: P+H* = H-P+, (2.21a) 

P 31 P 32 H a = H+P 3l P 32 , P+P+H+ = H^P+P^. (2.21b) 

That is, in Class (1, 1) we have considered a set of intermediate Hamiltonians H 11 and H sl 
both at the place between P^ and P 3 % and at the place between P 3 ^ and P 3 % simultaneously 
while in Class (0,1) and Class (1,0) we have considered an intermediate Hamiltonian at 
only one intermediate place, at the place between P^ and P 3 ^ in the former and at the place 
between P^ and P 3 % in the latter. 

But we easily find that the set of the intertwining relations (2.21) in Class (1,0) is 
transformed to the one (2.20) in Class (0, 1) by the transformations (2.18) which are caused 
by W(x) — > — W(x) provided that if 11 is transformed to H sl simultaneously under the same 
transformation. We shall hereafter call the set of the transformations (2.18) accompanied 
with the interchange H sl <H> H 11 the reflective transformation of a type A 3-fold SUSY 
system. Hence, as long as only one set of intermediate Hamiltonians is concerned, we only 
need to consider Class (0, 1) without any loss of generality. When we shall examine more 
than one sets of intermediate Hamiltonians, however, we must consider simultaneously both 
Class (0, 1) and Class (1, 0), as we shall discuss in Sections V and VI. 



III. EXISTENCE OF INTERMEDIATE HAMILTONIANS 

In this section, we shall derive the necessary and sufficient conditions for the existence 
of (at least) one set of intermediate Hamiltonians in each class classified in the last section. 
We note that Class (1, 1) can be regarded as a special case of either Class (0, 1) or Class 
(1,0). Thus, we only consider models in Class (0, 1) and Class (1, 0) which do not belong to 
Class (1,1) without any loss of generality. In the subsequent sections, we shall investigate 
each case in details. 



A. Conditions for Class (1,1) 

The necessary and sufficient condition for satisfying the first formula (2.19a) is that there 
exists a constant C 33 such that H and H l1 are expressed as 

2H~ = P 3 + 3 P 33 + 2C 33 = -d 2 -W - E' + W 2 + 2EW + E 2 + 2C 33 , (3.1a) 
2H [1 = P 33 P+ + 2C 33 = -d 2 + W + E' + W 2 + 2EW + E 2 + 2C 33 . (3.1b) 

Similarly, the necessary and sufficient condition for satisfying the second formula in (2.19b) 
is that there exists another constant C 32 such that H l1 and H ]1 are expressed as 

2H' a = P+P 32 + 2C 32 = -d 2 -W' + W 2 + 2C 32 , (3.2a) 
2H ]1 = P 32 P+ + 2C 32 = -8 2 + W' + W 2 + 2C 32 , (3.2b) 
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and the necessary and sufficient condition for satisfying the third formula in (2.19c) is that 



there exists another constant C 3 \ such that H il and H + are expressed as 

2// jl = P+P 31 + 2C 3 i = -d 2 -W + E' + W 2 - 2EW + E 2 + 2C 31 , (3.3a) 

2H + = P 3 - 1 P 3 + + 2C 3 i = -d 2 +W - E' + W 2 - 2EW + E 2 + 2C 31 . (3.3b) 

Comparing (2.1a), (3.1), (3.2), and (3.3) each other, we obtain 

6(W + EW) + E' + E 2 = —Q(R + C 33 ), (3.4a) 

2(W + EW) + E' + E 2 = 2{C 32 - C 33 ), (3.4b) 

2(W + EW) -E' -E 2 = 2(C 3 i - C32), (3.4c) 

6(W^' + EW^) - -E 2 = 6(C , 3i + ^)- (3.4d) 

Arranging the set of formulas (3.4) and using the relations (2.6), we have 

2Q'(z) = C 33 - C 31 , A"(z) = -3(2/2 + C 33 + C 3 i) = 2C 32 - C 33 - C 31 . (3.5) 



In terms of the coefficients of polynomials A(z) and Q(z) in (2.4), the necessary and sufficient 
conditions (3.5) are written as 

a 4 = a 3 = b 2 = 0, 2a 2 = 2C 32 - C 33 - C 31 , ^ g ^ 

26i = 6*33 — C31, —3/? = C33 + C32 + C31. 

The last three equalities in (3.6) just determine the parameters C31, C32, and C33 in terms 
of the parameters of the original type A 3-fold SUSY system as 

C 3 i = -f -&i-/?, C 32 = ^-R, C 33 = -^ + h-R. (3.7) 

We note that the first equality in (3.6) is identical to the solvability condition of type A M- 
fold SUSY systems (cf., Ref. [14], Eq. (6.13)), namely, the necessary and sufficient condition 
for the pair of any type A A/"-fold SUSY Hamiltonians H ± to be not only quasi-solvable 
but also solvable. Hence, a type A 3-fold SUSY system with a set of two intermediate 
Hamiltonians H l1 and H l2 is always solvable, and conversely a solvable type A 3-fold SUSY 
system always admits a set of two intermediate Hamiltonians. 



B. Conditions for Class (0, 1) 

To consider the first formula in (2.20a), we first note that the second-order operator 
-P32-P33 belongs to a type A 2-fold supercharge 

P32P33 = (d + W x - ^j(d + W x + |) , (3.8) 

with 

W 1 (x) = W(x) + -¥-. (3.9) 
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Hence, the necessary and sufficient conditions for satisfying the first formula (2.20a) are that 
there exists a second-degree polynomial Qi(z) in z such that (cf., Ref. [6]) 

Q^z) = -z\x)W l {x) ) = 0, (3.10) 

and that H~ and W 1 are expressed as 

2H = "d^ + w * ~ T + T ~ 2Rl ~ 2W *' (3 ' lla) 

2W l = --^ + W^ 2 - y + ^ - 2R l + 2W[, (3.11b) 

where Ri is a constant. The necessary and sufficient condition for satisfying the second 
formula (2.20b) is that there exists another constant C 3 i such that H ]l and H + are expressed 
as 

2W 1 = P+Pn + 2C 31 = -d 2 — W' + E' + W 2 — 2£W/ + £ 2 + 2C 31 , (3.12a) 
2#+ = P 3 -P 3 + + 2C 31 = -9 2 + W' — E' + W 2 — 2EW + E 2 + 2C 31 . (3.12b) 

Comparing (2.1a), (3.11), and (3.12) each other, and using (3.9), we obtain 

Q(W + EW) -E' -E 2 = -12(R - i?i) = 4(C 3 i + Ri) = 6(<7 3 i + (3.13) 

which is equivalent, in view of the relations (2.6), to 

-6Q'(z) - A"(z) = -12(R - R x ) = 4(C 3 i + = 6(C 3 i + R). (3.14) 

The condition (3.10), if combined with (2.5), leads to 

A'(z) = 2[Q(z)-Q 1 (z)], (3.15) 

which means that A'(z) is also a polynomial of at most second degree. In terms of the 
coefficients of polynomials A(z) and Q(z) in (2.4), the necessary and sufficient conditions 
(3.14) and (3.15) are written as 

a 4 = a 3 + 26 2 = 0, a 2 + 3&i = -2(C 31 + R ± ), 3R = 2R 1 -C 31 . (3.16) 

The last two equalities in (3.16) just determine the constants C 3 i and Ri in terms of the 
parameters of the original type A 3-fold SUSY system as 

C 31 = -- -h-R, R^-^-^ + R. (3.17) 
3 6 2 v ; 

We now find that a model in Class (0, 1) which satisfies (3.16) can also satisfy the condition 
(3.6) and thus belong to Class (1, 1) if and only if a 3 = b 2 = 0. Hence, we assume a 3 b 2 7^ 
without any loss of generality in the subsequent analyses of Class (0, 1). With the latter 
assumption, however, the solvability condition cannot be satisfied inevitably. Hence, all the 
models in Class (0, 1) are only quasi-solvable but may not be completely solvable. 

Finally, we note that the superHamiltonian and type A 2-fold supercharges constructed 
from H~, H ]1 , and P 32 P 33 form a type A 2-fold superalgebra since H~ and H sl are a type 
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A 2-fold SUSY pair with respect to the operator P^P-i?,- 111 particular, the anti-commutator 
of the type A 2-fold supercharges constructed in this way reads in components 

P&P&P*Pn = ASp l \H~ + R), P 32 P 33 P^ 2 = AS^iW 1 + P), (3.18) 
where S^ ' 1 ^ is a monic polynomial of second-degree given by 

S^\t) =t 2 -(^ + b 1 )t- 2( fll - 26 )6 2 - ^(a 2 - 3&0a 2 . (3.19) 

C. Conditions for Class (1,0) 

The analysis for Class (1,0) is almost the same as for Class (0, 1) in the previous section. 
That is, the first formula in (2.21) requires that H and H l1 must form an ordinary SUSY 
pair with respect to P 33 : 

2H~ = P£P 33 + 2C 33 = -d 2 -W - E' + W 2 + 2EW + E 2 + 2C 33 , (3.20a) 
2H 11 = P 33 P+ + 2C 33 = -d 2 +W + E' + W 2 + 2EW + E 2 + 2C 33 , (3.20b) 

while the second formula in (2.21) requires that H a and H + must form a type A 2-fold 
SUSY pair with respect to P^P^ and its conjugate: 

E' E 2 

2H 11 = -d 2 + W 2 H — 2P 2 - 2W^ (3.21a) 

E' E 2 

2H + = -d 2 + W 2 2 - — + — - 2R 2 + 2W^ (3.21b) 

where W 2 (x) = W(x) — E(x)/2. By following a similar route as in the previous section, it 
is straightforward to show that the above requirements are satisfied if and only if 

a 4 = a 3 -2b 2 = 0, C 33 = -j + h-R, R 2 = + ^ + R. (3.22) 

The last two equalities in the above just determine the constants C33 and R 2 and thus only 
the first formula stands essentially as the necessary and sufficient condition for a system 
to belong to Class (1,0). As in the case of Class (0,1), we assume a 3 b 2 7^ without any 
loss of generality to prevent a model in Class (1,0) from belonging also to Class (1,1). 
Hence, by the assumption, any system in Class (1,0) is only quasi-solvable but may not 
be completely solvable. In addition, the superHamiltonian and type A 2-fold supercharges 
constructed from iP 1 , H + , and P31P32 form another type A 2-fold superalgebra since H a 
and H + are a type A 2- fold SUSY pair with respect to the operator P 3 ~^P 3 ^. In particular, 
the anti-commutator of the type A 2-fold supercharges constructed in this way reads in 
components 

P&P&PfiPn = ^'°\H [1 + P), P 31 P 32 P 3 + 2 P 3 \ = ASf> l \H+ + P), (3.23) 
where S^' ^ is a monic polynomial of second-degree given by 

S£'°\t) =t 2 -(^-b i y + 2(ai + 26o)6 2 - l(a 2 + 3&i)a 2 . (3.24) 
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IV. CLASSIFICATION: CASES OF ONLY ONE SET 



We are now in a position to classify completely all the possible type A 3-fold SUSY 
potentials which admit one and only one set of intermediate Hamiltonians. In the previous 

section, we showed that there are three different classes, Class (1, 1), Class (0, 1), and Class 

(1,0), and that the necessary and sufficient conditions for the existence of intermediate 
Hamiltonians are 

Class (1, 1): a 4 = a 3 = b 2 = 0, (4.1) 

Class (0, 1): a 4 = a 3 + 2b 2 = (a 3 b 2 ^ 0), (4.2) 

Class (1, 0): a 4 = a 3 - 2b 2 = {a 3 b 2 ^ 0), (4.3) 

where we have neglected the other irrelevant constraints on the parameter relations. It is 
evident that these conditions explicitly break the GL(2, C) covariance of the original type A 
3-fold systems. Hence, we cannot apply the same classification scheme as the one in Ref. [6] 
which employs the full GL(2,C) covariance of the systems. However, we note that there 
exists a residual symmetry remained intact in the present cases. In fact, the transformation 
formulas (2.12) for the functions E(x) and W(x) tell us that for the set of linear projective 
transformations with 7 = the functions E(x) and W(x) are both invariant. The latter fact 
means in particular that each factor of the type A 3-fold supercharges (i = 1, 2, 3) in (2.2) 
and thus the intermediate Hamiltonians H a and ffl 1 introduced through the intertwining 
relations (2.19)-(2.21) as well are all invariant under an arbitrary GL(2, C) transformation 
with 7 = 0. The set of linear projective transformations with 7 = consists of the set of 
(complex) inhomogeneous linear transformations 

z = aw + /3 (a,^C,a/0), (4.4) 

where we have set 5 = 1 without any loss of generality. For an arbitrary inhomogeneous 
linear transformation (4.4), the transformation matrices of the parameters a>i {% — 0, . . . ,4) 
and bi (i = 0, 1,2) given by (2.10) and (2.11) respectively become triangle, and all of the 
conditions (4.1)-(4.3) are covariant, that is, 

«4 = a 3 = b 2 = =>- a 4 = a 3 = b 2 = 0. (4-5) 
a 4 = a 3 ±2b 2 = a 4 = 0, a 3 ± 2b 2 = a{a 3 ± 2b 2 ) = 0. (4.6) 

Therefore, all type A 3-fold SUSY systems with one set of intermediate Hamiltonians are 
classified by considering the equivalence class under the set of inhomogeneous linear trans- 
formations (4.4). We now easily show that the representatives of A(z) under this equivalence 
class in each case can be chosen, by noting the constraint a 3 7^ in Class (0, 1) and Class 
(1, 0), as listed in Table I. We note that Case II and Case IF (Case IV and Case IV, respec- 
tively) are transformed to each other by a projective transformation of GL(2, C) but not 
by any inhomogeneous linear transformation (4.4). That is the reason why we must treat 
them as separate cases though they are classified as one equivalent case in the classification 
of general type A A/"-fold SUSY models [6] . The still arbitrarily determinable constant a in 
Table I can be fixed by considering the scaling relations under the scale transformations of 
the parameter space spanned by {a^, b iy R}: 

z(x; vdi, ubi, uR) = z[ypux\ di, bi, R), (4.7a) 
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Case Class (1,1) Class (0, 1), (1, 0) 



I 


a/2 




II 


2z 




ir 




2z 3 


in 


z 2 /2 




IV 


2a{z 2 - 1) 




IV 




2az 2 (z+l) 


V 




2z 3 - g 2 z/2 - g 3 /2 



TABLE I: The complex classification scheme of type A 3-fold SUSY models with one set of in- 
termediate Hamiltonians. In the above, 02,53 £ C satisfy <7 2 3 7^ 27 g 2 and a G C is an arbitrary 
constant. 



E(x; vdi, ubi, uR) = y/uE(y/ux; a iy b iy R), (4-7b) 
W(x; vdi, vbi, vR) = \ZuW(y/ux; a i: b i: R), (4- 7c) 

V(x; vdi, ubi, uR) = vV{ypux\ a iy b iy R). (4-7d) 

Hence, in what follows we set a = 1 in all the cases without any loss of generality. For the 
classification, we recall the fact that Class (0, 1) and Class (1, 0) are connected with each 
other by the reflective transformation. Hence, we shall present only models belonging to 
Class (0, 1). We shall show in each case the change of variable z = z(x) and the two functions 
E(x) and W(x) determined by (2.5) as well as the potential parts of the 3-fold SUSY pair 
Hamiltonians and of the intermediate Hamiltonian(s) H l1 and/or H sl determined by 
(3.1)-(3.3) in Class (1, 1) and by (3.11)-(3.12) in Class (0, 1). 



A. Classification of Class (1,1) 

I) A{z) = 1/2. 
Functions: 



Potentials: 



z(x) = x, E(x) = 0, W(x) = -hx - 60. (4.8) 



2V~(x) = (b lX + b ) 2 + 3h - 2R, (4.9a) 

2V il (x) = (b lX + b ) 2 + b x -2R, (4.9b) 

2V ]1 (x) = (hx + bo) 2 -61 -2R, (4.9c) 

2V + (x) = (hx + bo) 2 -36i -2R. (4.9d) 



In this case, all the potentials (4.9) are harmonic oscillators. 
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II) A(z) = 2z. 
Functions: 



Potentials: 



z[x) = x , 



E(x) 



x 



W(x) = --x--. 



V-(x) = ^x 2 + 



hi 2 , (6o-2)(6q-4) , hibo + 3) 



8x 2 



+ 



R, 



6 (6 - 2) 6i(6o + l) 



+ 



V^(x) = -±x 2 + 
8 

8 



bl 2 6o(6o + 2) , 6l(6o-l) 



ix z 



+ 



R, 



R, 



(6 + 2)(6q + 4) 6i(6o-3) 
8x 2 4 



(4.10) 

(4.11a) 
(4.11b) 
(4.11c) 
(4.11d) 



In this case, all the potentials (4.11) are radial harmonic oscillators with centrifugal 
potentials. 



Ill) A(z) = z 2 /2. 
Functions: 



Potentials: 



z(x) 



E(x) 



W(x) = -61 - 6 e" 



2V~(x) = 6 (26i - 3)e~ x + 6 2 e 



2„-2x 



2R, 



2V a (x) = 6 (26i - l)e~ x + b 2 e~ 2x - 2R, 
2\/ jl (x) = 6 (26i + l)e~ x + b 2 e~ 2x - 2R, 
2V + (x) = 6 (26i + 3)e~ x + b^e~ 2x - 2R, 



(4.12) 



(4.13a) 
(4.13b) 
(4.13c) 
(4.13d) 



where R = R — 6 1 2 /2 — 1/3 is a constant. In this case, all the potentials (4.13) are Morse 
potentials. 



IV) A(z) = 2(z 2 - 1). 
Functions: 

z(x) = cosh 2x, 

Potentials: 

V~{x) 
V l \x) 
V' s \x) 



E(x) 



2 cosh 2x 
sinh2x 



W(x) = 



bi cosh 2x + 6 
2 sinh 2x 



26 (6i - 6) cosh 2x + 6 2 + (61 - 4) (61 - 8) 
8 sinh 2 2x 

26 (61 - 2) cosh 2x + 6 2 + 61 (6 X - 4) - 

2 > 

8 sinh 2x 

26 (6i + 2) cosh 2x + b 2 + 61 (61 + 4) - 
8 sinh 2x 
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R, 



(4.14) 

(4.15a) 
(4.15b) 
(4.15c) 



= 2W + 6)cosh2, + ^ + (5 1 + 4)(5 1 + 8) _ 
V ; 8sinh 2 2a; V ; 

where R = R — b 2 /8 — 4/3 is a constant. In this case, all the potentials (4.15) are Poschl- 
Teller potentials. We note that a system of Scarf potentials can be obtained by the choice 
A(z) = 2(z 2 + 1) which is connected with the representative A(z) = 2(z 2 — 1) by the 
combination of a complex linear transformation (4.4) and a scale transformation (4.7). 

It is interesting to observe that all the potentials given by (4.9), (4.11), (4.13), and (4.15) 
are shape invariant as can be easily checked by some obvious scaling of parameters [15]. It 
bears mention that the schemes of SUSY quantum mechanics [15-18] are consistent with the 
factorization method [19], intertwining relationships [20], and the shape invariance condi- 
tion [21]. The latter was first utilized by Gendenshtein [21] to track down solvable potentials. 
To be shape invariant the partner potentials while sharing a similar coordinate dependence 
can at most differ in the presence of some parameters as precisely has happened in the 
potentials above. 



B. Classification of Class (0, 1) 

II') A(z) = 2z 3 . 
Functions: 

z(x) = ^ e(x) = -1 w ^ = -Yx + h i x+h i x '- (416) 

Potentials: 

V~{x) = + -^x 4 + 1 V + — -b 1 -R, 4.17a 

8 4 8 8x z 

V 1 (x) = —x + —x + x +-^-R, (4.17b) 

V + (x) = -fx 6 + -^V + 1 g V + ^ + i~ R - ( 4 - 17c ) 

In this case, all the potentials (4.17) are well-known quasi-solvable sextic anharmonic 
oscillators [22]. 

IV) A(z) = 2z 2 {z + l). 
Functions: 

z(x) = ^-, E(x) = - 3+ l 2sinh * X , w{x) = ^* + bi*f*-\ (4 . 18) 
sinh x smhxcosha; 2 smh x cosh x 

Potentials: 

V~(x) = — cosh 4 x + - — - — cosh 2 x 

8 8 

_ (6l _ 6o + 3)(6 1 -6 + 5) + ^ + 6o_ fi| 
8 cosh x 8 sinh x 2 



13 



V* (*) = £ cosh 4 x + (26l - 3&0 + 4)b ° cosh 2 X 
V ; 8 8 



(6 1 -6 -l)(6i-6 + l) 



+ 



35 



2 ~~ 7T - ^> 



8 cosh^ x 8 sinh^ x 2 

K>) = ^ cosh', + < M '- 36 ° + 12 > 6 ° cosh' 
V ; 8 8 



X 



(6i-6 -3)(&i-6q-1) 
8 cosh 2 x 



+ 
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8 sinh x 



-b -R, 



(4.19b) 



(4.19c) 



where R — R — (b± — b )(bi — 36o)/8 — 4/3 is a constant. In this case, all the potentials 
(4.19) are quasi-solvable deformed Poschl-Teller or Scarf potentials [23]. 



V) A(z) = 2z 3 - g 2 z/2 - g 3 /2. 
Functions: 



z(x) = p(x), E(x) 



P"{x) 
p'(x) ' 



W(x) 



p{x) 2 - bipjx) - b 
p'jx) 



(4.20) 



Potentials: 
V-(x) = 

V }1 jx) = 

V + (x) = 



bipjx) + b 
2p'jx) 2 

hpjx) + b 
2p'jxf 

bipjx) + b 
2p'jxf 



b 1 p(x) + b - yp"(i) 
2 

hpjx) + b + -p"jx) 
g 

b 1 p(x) + b + -p"jx) 



+ ^p(x)-^p(2x)- b j-R, (4.21b) 
+ ™p(x) + ^p(2x)-^-R. (4.21c) 



In the above, pjx) is the Weierstrass elliptic function and bp = bp — g 2 /12 is another pa- 
rameter. The first term of each the potential is a rational function of p(x) by the formulas 
p'jx) 2 = 4p(x) 3 — g 2 pjx) — <?3 and p"(x) = Qpjx) 2 — g 2 /2. In this case, all the potentials 
(4.21) are quasi-solvable one-body elliptic BC type Inozemtsev potentials (cf., Ref. [24] and 
the references cited therein). 

V. DIFFERENT SETS OF INTERMEDIATE HAMILTON! ANS 



In this section, we shall study under what conditions a type A 3-fold SUSY system can 
have more than one sets of intermediate Hamiltonians. The latter possibility originates from 
the fact that each factor of a type A A/"-fold supercharge in a factorized form is not invariant 
under a subset of the GL(2, C) transformations although any type A A/"-fold supercharge 
as a whole is invariant under all the GL(2, C) transformations [6]. For the Af = 3 case, we 
easily check the latter fact directly from the definition (2.2) and the transformation formulas 
(2.12): 



2jz' 
72: — a 



P±[W,E] = P±[W,E] 



2 1Z ' 
•yz — a 



(5.1) 



Hence, the factors P 3 ^ and P 3 % in P^ are in fact not invariant under any GL(2, C) trans- 
formation with It means in particular that intermediate Hamiltonians H 11 and/or 
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H ]1 defined by the relations (2.19), (2.20), or (2.21) could not be invariant under any trans- 
formation with 7^0, even if they exist after the transformation. In fact, the existence 
of intermediate Hamiltonians after a transformation is not guaranteed automatically since 
the conditions for the existence (4.1) or (4.2) are not preserved under any GL(2,C) trans- 
formation with 7 7^ as was already discussed in Section IV. Therefore, another different 
set of intermediate Hamiltonians exists if there is a GL(2, C) transformation with 7^0 
for which the transformed parameters a« and hi given by (2.10) and (2.11) also satisfy the 
existence conditions (4.1) or (4.2). Furthermore, if there exist simultaneously such n(> 1) 
GL(2,C) transformations characterized by the sets {a i: 7*, Si}f =1 satisfying 7« 7^ and 
a i/li 7^ a j/lj f° r & U i 7^ h then the corresponding n sets of intermediate Hamiltonians are 
different from each other since the formula (2.12) tells us that the deformation is character- 
ized by the one parameter a/7. We shall hereafter say that two GL(2,C) transformations 
with 7j,7j 7^ are inequivalent if 0:1/7* 7^ a j/lj- 

In the classification of the systems which admit one set of intermediate Hamiltonians 
in Section IV, we considered the two different classes, namely, Class (1,1) where an in- 
termediate Hamiltonian exists at each of the two intermediate positions in the factorized 
type A 3-fold supercharge and Class (0, 1) where an intermediate Hamiltonian exists only 
at one intermediate position. And any system belonging to Class (1,0) can be obtained 
from the corresponding system belonging to Class (0, 1) by the reflective transformation 
W(x) —7- —W(x). Accordingly, we can consider a type A 3-fold SUSY system which admits 
simultaneously m different sets of intermediate Hamiltonians {H lk , ffl k }™ =1 of Class (1,1) 
and n different sets of an intermediate Hamiltonian {H ik }^^ +1 of Class (0, 1). We shall 
call such a class of systems Class (m,m + n) with an obvious implication of the terminology. 
It is evident that any system which admits simultaneously m different sets of intermediate 
Hamiltonians {H lk , Hi k }™ =1 of Class (1,1) and n different sets of an intermediate Hamil- 
tonian {H lk }™=™ +1 of Class (1,0), which would be called to belong Class (m + n,m), can 
be obtained from the corresponding system belonging to Class (m, m + n) by the reflective 
transformation (with the obvious accompaniment of the interchanges H^ k -H- H lk for k > 1). 
To investigate each class systematically, we first note that any system which belongs to Class 
(to, m + n) with m, n > always belongs to Class (1, 1) since the conditions for Class (1, 1) 
is stricter than for Class (0, 1). In other words, an arbitrary system in Class (to, m + n) with 
to > is a special case of one of the four systems in Class (1,1) classified in Section IV A 
and thus is always solvable in particular. Hence, only Class (0,n) with n > 1, which is an 
abbreviation for Class (0, + n), can have a quasi-solvable system which must be a special 
case of one of the three systems in Class (0, 1) classified in Section IV B. 

In addition to these classes, there could exist hybrid classes. That is, there is the possi- 
bility that a type A 3-fold SUSY system admits simultaneously different sets of intermediate 
Hamiltonians some of which belong to Class (0, 1) and the others of which belong to Class 
(1,0). These systems, if exist, could belong to neither Class (to, to + n) nor Class (0,n) 
since Class (0,1) and Class (1,0) are, as already mentioned, connected by the reflective 
transformation but not by any GL(2, C) transformation. We can thus consider a class of 
systems belonging to Class (to, m + n) which admit simultaneously I additional different 
sets of intermediate Hamiltonians {H lk Yk=L+i of Class (1,0). We shall call such a class of 
systems Class (m + l,m + n). We can assume n > I without any loss of generality since the 
reflective transformation interchanges Class (m + n, m + I) and Class (m + l,m + n). By 
following an argument similar to in the last paragraph, we shall separate the hybrid classes 
in two, the one is Class (m + l,m + n) with m,l,n > which is a special case of Class 
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(m,m + n) with m,n > and the other is Class (l;n), which is an abbreviation for Class 
(0 + I, + n), with n > I > 1 and is a special case of Class (0,n) with n > 1. 

In the subsequent sections, we shall study Class (m, m + n) with m, n > and Class 
(0, n) with n > 1 separately. 

A. Conditions for Class (m, m + n) with m, n > 

The necessary and sufficient conditions for the existence of a set of intermediate Hamil- 
tonians of Class (1,1) are given by (4.1). A system in Class (1,1) also belongs to 
Class (m, m + n) if there exist simultaneously n inequivalent GL(2, C) transformations for 
which the conditions (4.2) for Class (0, 1) are satisfied by the transformed coefficients a; 
(i — 0, . . . , 4) and 6, (i = 0, 1, 2), that is, 

a 4 = a 3 + 2b 2 = (a 3 6 2 ^0), (5.2) 

and in addition if there exist simultaneously m — 1 additional inequivalent GL(2, C) trans- 
formations for which the conditions (4.1) for Class (1, 1) are satisfied by hi and hi, that 
is, 

a 4 = a 3 = b 2 = 0. (5.3) 

From the transformation formulas (2.10) and (2.11), we see that when the condition (4.1) is 
satisfied, the transformed coefficients a 4 , a 3 , and b 2 respectively read 

A 2 a 4 = a 2 7 2 a 2 + a^ai + 7 4 ao, (5.4a) 
A 2 a 3 = 2cry(a5 + ^)a 2 + 7 2 (3a5 + £7)01 + 47 3 <5a , (5.4b) 
A6 2 = CK7&1 + l\- (5.4c) 

Hence, except for the trivial case 7 = 0, the conditions (5.2) are satisfied if and only if 

a 2 a 2 + a^ai + 7 2 a = 0, 2aa 2 + 701 — 2ab\ — 2^b Q = 0, (5.5) 

with a 3 b 2 7^ 0, and the conditions (5.3) are satisfied if and only if 

a 2 a 2 + cryai + 7 2 a = 0, 2aa 2 + 701 = 0, ab\ + 7& = 0. (5.6) 

In addition, a system of Class (m, m + ri) also belongs to Class (m + /, m + n) if there exist 
simultaneously additional / inequivalent GL(2, C) transformations for which the transformed 
coefficients hi and 6, satisfy the conditions (4.3) for Class (1,0), namely, 

a 4 = a 3 -26 2 = (a 3 6 2 ^ 0). (5.7) 

By the transformation formulas (5.4), they are satisfied if and only if (with a 3 b 2 7^ 0) 

a 2 a 2 + a^ai + 7 2 a = 0, 2aa 2 + 701 + 2ab\ + 2760 = 0. (5.8) 

In Section VIA, we shall investigate the conditions (5.5), (5.6), and (5.8) in each case 
separately. 



16 



B. Conditions for Class (0, n) with n > 1 



The necessary and sufficient conditions for the existence of an intermediate Hamiltonian 
of Class (0, 1) are given by (4.2). A system in Class (0, 1) also belongs to Class (0, n) if there 
exist simultaneously n — \ inequivalent GL(2, C) transformations for which the transformed 
coefficients di and bi satisfy the condition for Class (0,1), namely, Eq. (5.2). From the 
transformation formulas (2.10) and (2.11), we see that when the condition (4.2) is satisfied, 
the transformed coefficients a 4 , a 3 , and b 2 respectively read 

. O q OO Q A / \ 

A 04 = a 7<2 3 + a 7 a 2 + «7 ai + 7 ao, (5.9a) 
A 2 a 3 = a 2 (a5 + 3/?7)a 3 + 2cry(a<5 + ^)a 2 + 7 2 (3a<5 + £7)01 + 47 3 <5a , (5.9b) 
Ab 2 = -a 2 a 3 /2 + CK7&1 + >y\. (5.9c) 

Hence, we obtain 

A 2 (a 3 + 2b 2 ) = 4a 2 f3>ya 3 + 2cry(a<5 + ^)a 2 + 7 2 (3a<5 + £7)01 + 47 3 <5a 

+ 2Aa 7 &i + 2A 7 2 6 . (5.10) 

Therefore, the condition (5.2) is satisfied, except for the trivial case 7 = 0, if and only if 

a 3 a 3 + a 2 7a 2 + «7 2 ai + 7 3 a = 0, (5.11a) 
4a 2 a 3 + 2«7a 2 + 7 2 ai - 2cry6i - 27 2 6 = 0, (5.11b) 

with a 3 6 2 7^ 0. The second condition is derived by the elimination of ao- On the other hand, 
a system in Class (0,n) also belongs to Class (l;n) if there exist simultaneously I additional 
inequivalent GL(2, C) transformations for which the transformed coefficients di and bi satisfy 
the conditions for Class (1,0), namely, Eq. (5.7). By the transformation formulas (5.9) they 
are satisfied if and only if (with a 3 &27 7^ 0) 

a 3 a 3 + a 2 7a 2 + cry 2 ai + 7 3 a = 0, (5.12a) 
2a 2 a 3 + 2«7a 2 + 7V + 2a-fb x + 2-f\ = 0. (5.12b) 

In Section VI B, we shall investigate the conditions (5.11) and (5.12) in each case separately. 

VI. CLASSIFICATION: CASES OF MORE THAN ONE SETS 

Now that we have derived the existence conditions for another different set of intermediate 
Hamiltonians in a general form, we are in a position to proceed a detailed analysis for each 
case classified in Section IV. In what follows, we first investigate the systems which belong 
to Class (m,m+n) with m, n > and next the ones which belong to Class (0, n) with n > 1. 
All the former models are not only quasi-solvable but also solvable since they are special 
cases of Class (1,1). On the other hand, all the latter models are merely quasi-solvable since 
we have excluded either cases of a 3 = b 2 = or of a 3 = b 2 = in Class (0, n). 
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A. Classification of Class (m, m + n) with m, n > 



I) = 1/2: 

In this case, both the conditions (5.5) and (5.6) only have a trivial solution 7 = 0. Thus, 
the system admits no different sets of intermediate Hamiltonians. 

II) A{z) = 2z: 

In this case, the conditions (5.6) only have a trivial solution but the conditions (5.5) have 
one set of non-trivial solutions 

a = 0, 60 = 1 with a 3 = -2b 2 = 2 7 //3(^ 0). (6.1) 

Thus, the system admits another set of an intermediate Hamiltonian of Class (0, 1) and 
belongs to Class (1,2), which is an abbreviation for Class (1,1 + 1). On the other hand, the 
conditions (5.8) also have one set of non-trivial solutions 

a = 0, b = -1 with a 3 = 2b 2 = 0). (6.2) 

However, the latter GL(2,C) transformation a = is equivalent to the one in (6.1), and 
the latter solution b = — 1 is not compatible with bo = I in (6-1). Hence, the system does 
not admit a realization of Class (1 + 1,1 + 1). 



Functions: 



13 b 1 

z{x)=x 2 , E(x) = -, E{x) = — , W(x) = W(x) = ~x - — . (6.3) 

Ob Ob — £0b 

Potentials: 

V-(x) = b -tx 2 + JL + h-R, (6.4a) 

V-{x) = b p - ± 2 + I - R, V>\x) = I x 2 + A _ Ri (6 .4b) 

V*(x) = b p + ^-*, (6.4c) 

V+(x) = £a? + ^-^-R. (6.4d) 

It is interesting to note that the formulas (3.1) and (3.2) are not valid with E(x) and W(x) 
given in (6.3); the potential term V l2 (x), for instance, calculated by (3.1) and calculated by 
(3.2) do not coincide with each other. Hence, the system admits only H ]2 but not H l2 as 
the second set of intermediate Hamiltonians. All the potentials (6.4) including the newly 
appeared intermediate one V s2 {x) in (6.4c) belong to the class of radial harmonic oscillators 
with a particular angular momentum. 

Ill) A(z) = z 2 /2: 

In this case, only the conditions (5.6) have one set of non-trivial solutions 

a = 0, bo = with a 3 = 62 = 0. (6.5) 
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Thus, the system admits another set of intermediate Hamiltonians of Class (1,1) and belongs 
to Class (2,2). However, when the conditions (6.5) are satisfied, then E(x) = —E(x) = 1 
and W(x) = W{x) = —b\ so that all the potentials are just an identical constant. Hence, 
this case is nothing more than a trivial system. 

IV) A(z) = 2(z 2 - 1): 

In this case, the conditions (5.6) only have a trivial solution but the conditions (5.5) have 
two sets of non-trivial solutions 



a 



±7, b = Th±2 with a 3 = -2b 2 = -4 7 /(5 =f 0). 



(6.6) 



Thus, the system admits another set of intermediate Hamiltonians of Class (0, 1) correspond- 
ing to each of the solutions and belongs to Class (1,2). In addition, the system with the 
specific values of parameters b± — 2 and bo = admits the two different solutions simultane- 
ously and thus can have additional two different sets of intermediate Hamiltonians of Class 
(0, 1) corresponding to the two solutions. In the latter case, the system belongs to Class 
(1,3). On the other hand, the conditions (5.8) also have two sets of non-trivial solutions 



a = ± 7 , b = TbiT2 with a 3 = 2b 2 = -4 7 /(<5 =f p)& 0). 



(6.7) 



They are compatible with (6.6) if and only if b± — and b ± 2 with the inequivalent 
transformations, namely, a = ±7 for the former and a = =F 7 for the latter. In these 
particular two cases, the system belongs to Class (1 + 1, 1 + 1). 



IV-1) Class (1,2) 
Parameters: 

b = T(h-2), a 3 = -2b 2 
Functions: 



4 7 



Potentials: 



a 2 



2(0 ± 55) 



61 = -h + 



45 



5T(3 



z(x) = 
E(x) = 



V~(x 



V n (x 



V'\x 
V i2 (x 
V + (x 



cosh 2x, 

2 cosh 2x 
sinh2x ' 



W(x) = W(x) 



bi cosh 2x + b 







E(x) = 



2 cosh 2x 



2 sinh 2x 
4 sinh 2x 



sinh 2x cosh 2x =F 1 



T(6i - 2) (61 - 6) cosh 2x + b\ - 8&1 + 18 



4 sinh 2x 
T(h - 2f cosh 2x + b\ - 46i + 2 

4 sinh 2 2x 
=F(&! 2 -4) cosh2a; + 6 1 2 + 2 



-R, 



R, 



-R, 



8 



4 sinh 2x 
^f(bf -4) cosh2x + 6 1 2 + 2 

4 sinh 2 2x cosh 2x =F 1 

=F(6i - 2) (61 + 6) cosh 2x + 6 2 + 46i + 18 
4 sinh 2 2x 



-R, 
-R, 



(6.8) 



(6.9) 



(6.10a) 
(6.10b) 
(6.10c) 
(6.10d) 
(6.10e) 
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where R = R — b 2 /8 — 4/3 is a constant. The formulas (3.1) and (3.2) are again not valid 
with E(x) and W(x) given in (6.9) and thus the system admits only H ]2 but not H l2 as the 
second set of intermediate Hamiltonians. All the potentials (6.10) but the newly appeared 
intermediate one V l2 (x) in (6.10d) belong to the class of Poschl-Teller potential with only 
one free parameter. The deformation term in (6.10d) is reminiscent of the generalized 
Poschl-Teller potentials constructed in Ref. [4]. It is worth noting that the deformed 
potential (6.10d) is connected with the other shape- invariant potentials by the intertwining 
relations (2.20) and thus is almost isospectral to them and is in particular solvable. 



IV-2) Class (1, 3) 
Parameters: 



h = 2, a 3 = -2b 2 = 
b = 0, a 3 = -2b 2 = 



4 7 



5-/3' 
47 



Functions: 

z(x) = cosh2x, E(x) = 

E(x) 
Potentials: 



2 cosh 2x 
sinh 2x 



5 + p' 

2 cosh 2x 

sinh 2a; 
4 sinh 2x 



a 2 
a 2 



2Qg + 55) 
(3-5 ' 

2((3 - 55) 
(3 + 5 ' 



bi 

k 



W{x) = W{x) = W{x) = 
2 cosh 2x 



2(5 + (3) 
6-p ' 

2(5-13) 
5 + [3 ' 

cosh 2x 



(6.11) 



V~(x 



V n (x 



V' s2 (x 
V i3 (x 
V + (x 



cosh2x — 1 
3 

2 sinh 2x 



E(x) 



sinh 2x ' 
4 sinh 2x 



sinh 2x cosh 2x + 1 



-R, 



2 sinh 2 2x 
3 

+ 



R, V^(x) = 



2 sinh 2 2x 



R, 



8 



2 sinh 2 2x 
3 

2 sinh 2 2x 
15 

2 sinh 2 2x 



cosh 2x — 1 
8 

cosh2x + 1 



-R, 
-R, 



-R, 



(6.12) 

(6.13a) 
(6.13b) 
(6.13c) 
(6.13d) 
(6.13e) 



where R — R— 11/6 is a constant. There are essentially no free parameters in this case. 
1,1 + 1) 



IV-3) Class (1 
Parameters: 

61 = 0, a 3 = -2b 2 = 
b = ±2, h 3 = 2b 2 = 

Functions: 

z(x) = cosh2x, E(x) 

E(x) 



47 



5TP' 
47 



6±p' 
2 cosh 2x 



a 2 = 



a 2 



2((3± 55) 

PT5 
2(Pt$5) 



sinh 2x ' 
2 cosh 2x 4 sinh 2x 



P±5 

W(x) = W(x) = W(x 



61 = 
61 = 



sinh 2x cosh 2x =F 1 



E(x) 
20 



2 cosh 2x 
sinh 2x 



A5 
5^P' 
45 
= 5±f3' 

1 

~^ sinh 2x ' 
4 sinh 2x 



(6.14) 



(6.15) 



cosh 2x ± 1 



Potentials: 

V~(x 
V l \x 
V i2 (x 
V' i2 {x 
V + (x 



=F6 cosh 2x + 9 

2 sinh 2 2x 
=F2 cosh 2x + 1 



cosh2x ± 1 



2 sinh 2 2x 
=F2 cosh 2x + 1 

2 sinh 2 2x 
±2 cosh 2x + 1 

2 sinh 2 2x cosh 2x =f 1 
±6 cosh 2a; + 9 - 

2 > 

2 sinh 2x 



V ] \x) = 
-R, 



±2 cosh 2x + 1 
2 sinh 2 2x 



8 



R, 



(6.16a) 
(6.16b) 
(6.16c) 
(6.16d) 
(6.16e) 



where R = R — 4/3 is a constant. As in the previous case, there are essentially no free 
parameters in this case. 



B. Classification of Class (0, n) with n > 1 

II') A(z) = 2z 3 : 

In this case, the conditions (5.11) have one set of non-trivial solutions a = b = but with 
a 3 = b 2 = which should be excluded. In fact, the system corresponding to the latter 
solutions is identical with the one in the case II, which belongs to Class (1,2), already 
found in the previous section, Eqs. (6.4). 

IV) A{z) = 2z 2 {z + l): 

In this case, the conditions (5.11) have two sets of non-trivial solutions, the one is 

a = - 7 , b 1 = b -2 with a 3 = -26 2 = 2 7 /(/3 + (5)(^0), (6.17) 
but the other is 

a = b = with a 3 = b 2 = 0, (6.18) 

and thus should be discarded. Indeed, the choice of the latter solution simply leads to the 
system of Class (1,2), and together with the former solution (6.17), to the system of Class 
(1,3) which are identical with the systems (6.10) and (6.13) respectively in the case IV 
already found in the previous section. Hence, only the solution (6.17) leads to a new system 
which belongs to Class (0,2). 

On the other hand, the conditions (5.12) also have two sets of non-trivial solutions, the 
one is identical with (6.18) to be discarded while the other is 

a = -7, b 1 = b with a 3 = 26 2 = 2 7 /(/3 + 5)(^0). (6.19) 

Hence, the Class (0,1) system (4.19) admits an intermediate Hamiltonian of Class 
(1,0) and thus belongs to Class (1; 1). However, the GL(2,C) transformation a = —7 
of the latter solution is equivalent to the one in (6.17). Hence, the system does not 
admit a realization of Class (1;2). We note that the choice of the two solutions (6.18) 
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and (6.19) leads to the system of Class (1 + 1, 1 + 1) which is identical with the system 

IV'-l) Class (0, 2) 
Parameters: 

n - n t 2 7 - 2(5-2/3) t , 25 

b 1 = b - 2, a 3 = -2b 2 = ■ z , a 2 = — - ■ - , b x = -b 



(3 + 5 1 * 5 + (3 ' 1 u 13 + 5 
Functions: 

, . 1 . . ~ . . 6 sinh 2 x cosh 2 x — 2 sinh 2 x — 1 

z(x) = 5—, WMx) = WMx) = , 

sinhV 2sinhxcoshx 

. , 3 + 2sinh 2 x , 1 — 2 sinh 2 x 

E ^ = — T iT"' E ^ = iT"- 

sinh x cosh x smh x cosh x 

Potentials: 

V~ (x) = — cosh 4 x — — — — cosh 2 x ^-s 1 ^ 1 \~b Q — R, 

V ; 8 8 8 cosh 2 x 8 sinh 2 x 6 

V l\ x ) = b A cosh * x - b l cosh * x 3 + 35 + 11 _ ^ 

8 8 8 cosh x 8 sinh x 6 

V^s) = h l cosh 4 x - M cos h 2 x + ^r- + IT - ^ 

8 8 8 cosh x 8 sinh x 6 



M cosh 4 x - (&0 - cosh 2 x - 15 + 15 + 11 - h 
8 8 8cosh 2 x 8sinh 2 x 6 2 



IV-2) Class (1; 1) 
Parameters: 



0i = 0o, a 3 = 2b 2 = a , - , a 2 = — f , Q , »i = -0o 



Functions: 



(3 + 5' z ' 5 + (3 ' 1 u /3 + <T 



1 j 6 smh xcosh x-1 

z(x) = WMx) = WMx) = — — . 

sinh 2 x 2 sinh x cosh x 

. 3 + 2sinh 2 x - . 1 — 2 sinh 2 x 

E x = rr r— , £ x) = — — . 

smh x cosh x smh x cosh x 



Potentials: 






V~{x) = 


b 2 

— cosh 4 x — 
8 


(oo + 12)6, 
8 


V ]1 (x) = 


— cosh 4 x — 
8 


(oo - 4)6 
8 


V' l \x) = 


— cosh 4 x — 
8 


(&o + 4) oo 
8 


V + (x) = 


6 2 

— cosh 4 x — 




(oo - 12)6, 

c 



15 3 4 n 

" x 2— H 2~ + o + 77 - R > 

8 cosh x 8 sinh x 3 2 



,2 1 35 4 b n 

cosh x H 2 1 5 1 

8 cosh x 8 sinh x 3 2 



35 I 4 

8 cosh x 8 sinh x 3 
15 4 

8 cosh 2 x ' 8 sinh 2 x 3 



■ cosh 2 x — 1 : — ^ \~ - ~ b — R. 
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V) A(z) = 2z 3 -g 2 z/2-g 3 /2: 

In this case, the conditions (5.11) have three sets of non-trivial solutions 

a = e a , 6 = -e i 6 1 + 4e i 2 -^ with a 3 = -2b 2 = 0), (6.26) 

4 p - dd 

where each = p{oo-i) (i = 1, 2, 3) is the value of the Weierstrass elliptic function at the half 
of the fundamental period 2uji which satisfies the third-degree algebraic equation 

4e?-g 2 e t -g 3 = 0. (6.27) 

Thus, the system admits another set of intermediate Hamiltonians of Class (0, 1) correspond- 
ing to each of the solutions and in those cases the system belongs to Class (0, 2). In addition, 
the system with the specific values of parameters b\ = 4(ej + ej) and b = — (ej + ej) 2 — Se^ej 
(i ^ j) admits the two different solutions simultaneously and thus can have additional two 
different sets of intermediate Hamiltonians of Class (0, 1) corresponding to the two solu- 
tions. In the latter case, the system belongs to Class (0,3). We note, however, that the 
three different solutions are incompatible simultaneously and hence any Class (0,n) with 
n > 3 cannot be realized. On the other hand, the conditions (5.12) also have three sets of 
non-trivial solutions 

a = e tl , b = -e t b 1 -2e 2 + 9 j with a 3 = 2b 2 = 0). (6.28) 

Hence, the Class (0, 1) system (4.21) also admits an intermediate Hamiltonian of Class 
(1, 0) and thus belongs to Class (1; 1). In addition, a choice of one solution a = e^ in (6.26) 
and another a = e^ (j ^ i) in (6.28) is compatible with b\ = 2ej and b = e 2 — e^ej — e 2 . 
In the latter case, the system possesses two additional intermediate Hamiltonians the one 
belongs to Class (0, 1) and the other to Class (1,0), and thus it is a member of Class (1, 2). 
However, a choice of three solutions, e.g., two a = ej7, e/y in (6.26) and one a = e^ in 
(6.28) with i ^ j ^ k ^ i, conflict with the assumption of non-degeneracy g 2 3 ^ 27g 2 (for 
the latter example, = — 2ej must hold). Therefore, the hybrid classes such as Class (1,3) 
and Class (2; 2) cannot be realized anymore. 

V-l) Class (0,2) 
Parameters: 

h h +a 2 92 % , 4e?P-(12e? + g 3 )5 

b ° = ^ + ^ " 4 ' h = ~ h + 2e^-e t 5) ' 2q ^ 

i P><b . Ae 2 /3 2 + g 3 /3S-(Aef + g 3 ) ei S 2 { - ] 

a 3 — —zo 2 — — -, a 2 — 6- 



Functions: 



(3 -erf' 2e l ((3-e t 5) 2 

p(x) 2 - hp(x) - b 



z(x) = p(x), W(x) = W(x) = 



p'{x) 



w( v P"{x) . p"{x) 2p'{x) 

E(x) = ,, : , E(x) 



(6.30) 



p'(x) ' ' p'(x) p(x)-ei 
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Potentials: 
V~{x) = 
V ]1 (x) = 
V' i2 {x) = 



hpjx) +b 

2p'(xf 
bipjx) + bp 

hp{x) +b 
2p\xf 
2p> i ) 



+ 



hp(x) +b - yp"(l) 

2 

hp(x) + b + -p"(x) 
hp(x) +b + -p (x) 



-lp(x) + ^p(2x) + ^-R, (6.31a) 



40 



5 



+ ^ p{x) ~ i8 p(2a;) ~i- R > (6 - 31b) 



4 , . 5 . . 



V + (x) = 



p(x 
hpjx) + b 
2p\xf 



+ 4 ei -^-R, 

6i 2 



(6.31c) 



hp(x) +b + -p"(x) 



10 , . 55 . 36i „ 

+ irfpfr) + t^p( 2x ) ~ - R> ( 6 - 31d ) 



9 



18 c 



where bo = b — gij\2 = — e^i + 2p"(oji)/3. The first term of each the potential can be 
expressed solely in terms of p(x). With the latter value of bo, we have 



h^±^ [bMx) + - bo + C p'' { x)} 

[3b 1 (p(x)-e l ) + 2p"(u i )] 2 Cy, 3h(p(x) - e t ) + 2p"(co t ) 
72UUp(x) - ei ) 12^ 

where and hereafter i ^ j ^ k ^ i k — 1, 2, 3). 



p(x) - e { 



(6.32) 



V-2) Class (0, 3) 

Parameters: hi and hi are the same as (6.29). 

bi = -Ae k , % = -el - 3^-, bi = ~h + 



Functions: 



E{x) 
Potentials: 
V~{x) -- 
V ] \x) -- 



4eff3-(12ef + g 3 )5 
2e 3 {(5-e ] 5) ' 



a 3 = -2b 2 



(5 - e 3 8 ' 



p'(x) 

bipjx) + b 
2p'jxf 

hpjx) + bp 
2p\xf 



a 2 = 3 '■' 



4e//3 2 + g 3 /35 - (4e| + g,)e^ 



z(x) = p(x), W(x) = Wjx) = Wjx) = 

2p'(x) 



2e j (/3-e j 8)* 



p{x) 2 + 4e k pjx) + el - 3^ 
p'jx) 

p"jx) 2p'jx) 



t(x) p(x) - e. 



10 



, E{x) 



j'jx) p(x) - e 3 



hpjx) +b - —p"jx) 
2 

bip(x) +b + -p"jx) 



(6.33) 



(6.34) 



" lp(x) + ^p(2x) + ^-R, (6.35a) 

40 ^ h 

^ -P(x) - -p(2x) - j - R, (6.35b) 
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V ]2 jx) 



hpjx) + b 
2p>i) 



hpjx) + b + -p"jx) 
h 



p(x) - 6i 2 



V i3 jx) = 



hpjx) + b 
2p'jxf 



+ 



p(x) 



V + jx) 



hpjx) + bp 
2p'jxf 



hpjx) +b + -p"(x) 
+ 4e ? - - — -R, 



hpjx) + b + -p"jx) 



(6.35c) 



(6.35d) 

in oi, 

+ yP(*) + ^p(2x) - ^ - i2. (6.35e) 



With the values of &i and bp in (6.33), the first term of each the potential can be expressed 
solely in terms of p(x) as 



bipjx) + bp 
2p'jxf 



[hpjx) + b + Cp"jx)] 



2[3e k p(x) + el + 2eiej] 2 C ^ 3e k pjx) + e| + 2e^^ 

i=i 



V-3) Class (1; 1) 
Parameters: 





bp = -eibi 




a 3 = 2b 2 = 


Functions: 


z(x) = 
E(x) = 


Potentials: 




V'jx) = 


hp(x) + bp 
2p'jxf 


V ll jx) = 


bipjx) + bp 
2p'jxf [ 
2p"ju t ) 
pjx) - ei 


V ]1 jx) = 


bipjx) +b 
2p'jxf 


V + (x) = 


hpjx) + bp 
2p'jxf 



,2 92 
■i I - 4 , 



p"juj)l 
/9 - e,5 ' 



6i = -6i + 
a 2 = 3 



4e 2 /3 + (4ef + g 3 )<? 
2ei(/3 - e^) 
4ef/3 2 + g 3 p6 - j4 e ? + g 3 )e t 5 2 



p"jx) 
p'jx) ' 



Ejx) 



10 



2e,(/3 - e t 5) 2 

_ pjx) 2 - hpjx) - bp 

x ) ~ ,t \ ■, 

p'jx) 

p"jx) 2p'jx) 



Jjx) pjx)-e i ' 



hpjx) +b - —p"jx) 
4 

bipjx) +b - -p"jx) 
2 

hpjx) + b + -p"jx) 
g 

hpjx) +b + -p'jx) 



(6.36) 



(6.37) 



(6.38) 



8 Q1 3/)i 

- -Pjx) + -pj2x) + ^ - R, (6.39a) 

2 7 

" gp( x ) + ^P(2a:) 











(6.39b) 


40 . . 


— pj2x) - 

lg 5 V ) 


6i 
2 




(6.39c) 


10 . . 

+ -gPjx) + 


55 

18 p(2;r) " 


36i 
IT 




(6.39d) 
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where 60 = bp _ #2/12 = — e^i — p"((Ji)/3. The first term of each the potential can be 
expressed solely in terms of p(x). With the latter value of 60, we have 



2p'(xf 



[b 1 p(x) + b + Cp"(x)} 



12 ^ 

1=1 



72lK =1 (p(x)-e,) 
V-4) Class (1;2) 

Parameters: and fej are the same as (6.29). 



bi = 2ei, b = ef + eje k , b\ = —bi + 



p(x) - ei 



4e 2 (3 + (4ef + g 3 )5 
2e j ((5-e j 5) 



a 3 = 2b 2 



Functions: 

z{x) = p(x), W(x) = W(x) = W(x) 



p"(u,jh s _ Aef(3 2 + g 3 (35-(Aef + g 3 )e j 5 2 
p- ej 6' ° 2 2e j (p-e j 6) 2 

p(xf - 2e, L p{x) - e 2 - e 5 e k 



p'{x) 



E(x) ■■ 
Potentials: 
V~(x) = 

V l \x) = 



p"{x) 



, E(x) 



p"{x) 2p'{x) 



p'(x) 

bjp{x) + b Q 
2p\xf 

b x p(x) + b 
2p'{xf 



b lp (x) + 6 



+ 



V' sl (x) 
V' i2 {x) 



bipjx) + bp 
^ 2p'(xf 
bjpjx) + bp 
2p'{xf 
2p"(u t ) 



p'(x) p(x) - e 
4 „ 

b lp {x) + bp - -p (x) 
h 

+ 4ej + --R, 
2 

hp(x) +b + -p"(x) 
2 

hp(x) +b + -p (x) 



p"(x) 2p'(x) 



p'(x) p(x) - ej 



9 



7 



bi 



p(x) -6i 2 



V + (x) 



hpjx) + bp 
' WW 2 



b lP (x) +b + -p"(x) 



40 , , 5 /n N 
+ qP( x ) ~ -^P(2x) 



10 , v 55 /0 , 



36i 



(6.40) 



(6.41) 



(6.42) 



|p(a;) + ^p(2:r) + ^ - i2, (6.43a) 



(6.43b) 
(6.43c) 



(6.43d) 



- R. (6.43e) 



With the values of 61 and 6 in (6.41), the first term of each the potential can be expressed 
solely in terms of p(x) as 

hp{x) +b 



2p'{xf 



[b lP (x)+b + Cp"(x)} 
[3ejp(x) + e 2 + 2eje k ] 2 C_ y-> Sejp(x) + e 2 + 2e j e k 



(6.44) 
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VII. PARASUPERSYMMETRY AND GENERALIZED SUPERALGEBRAS 



In the case of J\f = 2, it was shown [1] that any type A 2-fold SUSY system which 
has an intermediate Hamiltonian admits a realization of second-order paraSUSY [9] and 
a generalized 2-fold superalgebra [10]. Thus, it is natural to ask whether an analogous 
realization is possible in the present M = 3 case. In what follows, we show that it is indeed 
the case. More precisely, a type A 3-fold SUSY system of Class (1,1) admits a realization of 
third-order paraSUSY [11, 12] and a generalized 3-fold superalgebra found in Ref. [13] while 
one of Class (0, 1) does only a restricted version of the latter algebra. We shall first discuss 
the former realization and then the latter. 

A. Parasupersymmetry of Order 3 in Class (1, 1) 

Higher-order paraSUSY was introduced in Refs. [11, 12] as a generalization of second- 
order one [9]. In the case of third-order, it is characterized by the following algebraic 
relations: 

(Q±) 3 ^ 0, (Q±) 4 = 0, [Q± H P ] = 0, (7.1a) 

(Qp) 3 <2p + (Qp) 2 QpQp + QpQp(Qp) 2 + QUQp) 3 = Q(Qp) 2 h p . (7.1b) 

By the introduction of parafermionic coordinates ijjp of order 3 satisfying [13] 

(#) 4 = 0, {^p,^} + {(^) 3 ,(#) 3 }=2/, {(^-) 2 ,(^) 2 }=J, (7.2) 

a quantum mechanical realization of paraSUSY of order 3 is achieved by defining the triple 
(H P , Q±) as 

Hp = Ho^-fdj+f + ifx(^p - (#) 2 (^p) 2 ) 

+ F 2 ((#) 2 (^p) 2 - (#) 3 (VV) 3 ) + # 3 (#) 3 (^p) 3 , (7.3a) 

Q P = Qi^pf^t? + Q~2 (V£(Vv) 2 - (#) 2 (^) 3 ) + Qs (#) 2 (^p) 3 , (7.3b) 

Q+ = Q+(^p) 2 (#) 3 + Qti^pi^tf - (^p) 2 (#) 3 ) + Q 3 + (#) 3 (^p) 2 , (7.3c) 

where 

Hk = ~l^ + Vk( > x) ' Qt = ^ + W k (x). (7.4) 

The linear space in which the system (Hp, Qp) shall be considered is J x V3 where 5 is a 
linear space of complex functions such as L 2 (R) and V 3 = X)fc=o^3^ * s ^ ne parafermionic 
Fock space of order 3 composed of the A;-parafermionic subspaces (k — 0, ... ,3). Then, 
the latter system satisfies the third-order paraSUSY algebra (7.1) if and only if [12, 13, 25] 

2H = Q-Q+ - 2R X , 2H X = Q+Q~ - 2R X = Q 2 Q\ - 2R 2 , (7.5a) 
2H 2 = Q\Q 2 - 2R 2 = - 2R 3 , 2H 3 = Q+Q7 - 2R 3 , (7.5b) 

where R^ (k = 1,2,3) are constants satisfying 

R 1 + R 2 + R 3 = 0. (7.6) 
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Comparing now the paraSUSY conditions (7.5) with the Class (1, 1) conditions (3.1)-(3.3), 
we immediately notice that any type A 3-fold SUSY system which belongs to Class (1, 1) 
admits a realization of paraSUSY of order 3 by the following identifications: 

Qt = P?4-k, Rk = -C 34 . k , H = H~, H l = H il , H 2 = H>\ H 3 = H + . (7.7) 

By the formulas (3.7) and the second relation in the above, the constraint (7.6) is expressed 
in terms of the type A A/"-fold SUSY parameters as R = 0, which is identical to the constraint 
in the M = 2 case (cf., Ref. [1], Section 4). 

The realization of third-order paraSUSY via the formulas (7.5) admits another nonlinear 
relation, Ref. [13], Eq. (4.33). For the present system, it reads by the formulas (3.7) and the 
second relation in (7.7) 

raw) 3 + ( SiHi^- 2 ) + 1 SiSri-H ) + (Qp) 3 (q 



k vQp~) 2 (QpTQp J I (QpT(Qp) 3 Qp j 

= 8^H P + R + ^) 2 -b^H P + R- 2 -fy (7.8) 

We note that this algebraic relation holds irrespective of the paraSUSY constraint (7.6). 
Thus, in this sense the latter algebra (7.8) is more general than the paraSUSY algebra 
(7-lb). 

In the subsector with the parafermion number and 3, the nonlinear algebra (7.8) reduces 

to 



{(Qrf, (Qp-) 3 } = 8 ((Hp + R + |) 2 - tf) (h p + R - ^ 



• (7-9) 

Sx(vf+Vf) 



This, together with the relations 

[(Q±f,H P ] = {(Q±) 3 ,(Q±) 3 } = 0, (7.10) 

which follow directly from the third-order paraSUSY relations in (7.1a) forms a 3-fold su- 
peralgebra. We can now easily check that the latter algebra exactly coincides with type A 
3-fold superalgebra (2.15) in Class (1, 1) by the conditions (4.1). Hence, an arbitrary type 
A 3-fold SUSY system which belongs to Class (1,1) admits a realization of paraSUSY of 
order 3 and the generalized type A 3-fold superalgebra (7.8). 



B. Generalized 3-fold Superalgebra in Class (0,1) and Class (1,0) 

Contrary to the case of Class (1,1), any system belonging to Class (0,1) and (1,0) 
admits neither paraSUSY of order 3 nor quasi-paraSUSY of order (3, q) [13]. The reason 
is the lack of a 'shape-invariant' condition at the place where an intermediate Hamiltonian 
is absent. However, as we shall show shortly, a restricted version of the generalized type A 
3-fold superalgebra (7.8) still holds in each of Class (0,1) and Class (1,0) with the same 
parafermionic setting as (7.3) and (7.7). Indeed, substituting the relations (7.7) into the 
formulas (4.21)-(4.26) in Ref. [13] and using the intertwining relations (2.20) and the formula 
(3.18) for Class (0,1), and (2.21) and (3.23) for Class (1,0), respectively, we see that in 
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the case of Class (0, 1) the following algebraic relation holds in the subsector with the 
parafermion number 0, 2, and 3 



(Q P ) 3 (Qp-) 3 + 



Q p (Q+) 3 (Qp) 2 
(QtnQ P ) 3 Qt 



H P + R) ( Hp + R + y + 6i 



Jx(vi ' +vf + V ( 3 >) 



(7.11) 



and that in the case of Class (1, 0) the algebra which holds in the subsector with the 
parafermion number 0,1, and 3 reads 



+^3 



(Qp) (Qp) 



Q+(Q P ) 3 (Q+) 2 
(Q p ) 2 (Q+) 3 Qp 



H I: >+ mlHr + R+j-h 



yx(vW+vW+vW) 



(7.12) 



In the whole parafermionic vector space J x V 3 , however, no algebraic relations like (7.8) 
hold for Class (0, 1) and Class (1, 0). In fact, for Class (0, 1) each of the second term in the 
l.h.s. of (7.8) is calculated in the subsector with the parafermion number 1 as 



Q+(Q P ) 3 (Q+) 2 = 2P-p+P+(Hi 1 - C 31 )P 3 2Uv«> 
(Qp) 2 (Qp1 3 Qp = 2P 3 + 2 (# jl - c 3 i)p 3 - 2 p 33 k 



(7.13) 



33| 5 xV 



and cannot be expressed as a polynomial of W 1 . Similarly, for Class (1, 0) each of the third 
term in the l.h.s. of (7.8) is calculated in the subsector with the parafermion number 2 as 



2 - f 31- f 31- P 32 {H l1 - C' 33 )P 3 t| JxV (2), 



Qp(Q+) 3 (Q p ) 2 
(Qp) 2 (Q p ) 3 Qp = 2P 32 (H l1 — C^P^P^P^ 



(7.14) 



and again cannot be expressed as a polynomial of H a . 

We note that both the nonlinear algebras (7.11) and (7.12) are compatible with the type 
A 3-fold superalgebra (2.15). In fact, it is easy to check that both (7.11) and (7.12) reduce 
to the anti-commutator of the type A 3-fold superalgebra (2.15b) in the more restricted 
subsector with the parafermion number and 3 by noting the condition (4.2) for the former 
Class (0, 1) and by (4.3) for the latter Class (1,0). 



VIII. DISCUSSION AND SUMMARY 

In this article, we have fully investigated the necessary and sufficient conditions for a type 
A 3-fold SUSY system to have one or more sets of intermediate Hamiltonians and then made 
the complete classification of them by the property of the GL(2, C) transformations. When 
only one set of intermediate Hamiltonians is concerned, there are three different patterns 
in the existence and called Class (1, 1), Class (0, 1), and Class (1,0), respectively. We have 
found that all the models which belong to Class (1,1) are not only solvable but also shape 
invariant while the ones which belong to Class (0, 1) or Class (1,0) are just quasi-solvable. 
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Case Possible Classes 



I 


(1,1) 






II 


(1,1) 




1,2) 


IF 


(0,1) 






III 


(1,1) 






IV 


(1,1) 




1,2) D (1 ' 3) 

1(1 + 1,1 + 1) 


IV 


(0,1) 




f(0,2) 

L(i; i) 


V 


(0,1) 


H 


f(0,2) D (0,3) 

[(1;1)d(1;2) 



TABLE II: The possible classes of intermediate Hamiltonians for each case of type A 3-fold SUSY 
models. 



When more than one sets of intermediate Hamiltonians are concerned, there emerge various 
patterns depending on the functional type of each model. In Table II, we summarize the 
possible classes for each case of type A 3-fold SUSY models. 

It is now evident from Table II that the structure of higher-order intertwining operators 
is much richer than the degree that one can classify them solely by the notion of reducibility 
introduced in Refs. [7, 8]. It is not only because the requirement of reality is restrictive 
but also because there are various patterns in the existence of intermediate Hamiltonians. 
Needless to say, the number of possible patterns drastically increase as the order M of 
intertwining operators gets higher. 

Although we have not assumed in this article the reality of Hamiltonians and thus have 
analyzed general complex Hamiltonians by employing the GL(2, C) transformations, it is 
straightforward to examine and classify real Hamiltonians by the use of the real GL(2, M) 
transformations instead of GL(2,C). We only show in Table III an example of the real 
classification scheme for that purpose. Note, however, that some of the possible classes for 
Case IV and Case V in Table II which can exist in the complex case might be missing in 
the real case since the solutions to the conditions (5.5), (5.6), (5.8), (5.11) or (5.12) are not 
necessarily real. 

The realization of the variant generalized 3-fold superalgebras in Section VII indicates 
that the parafermionic formulation like (7.3) could provide a more adequate and advanta- 
geous framework to formulate Af-fold SUSY than the conventional fermionic formulation 
like (2.14). In the conventional approach, the type A 3-fold superalgebra (2.15) cannot 
characterize nor detect the existence of intermediate Hamiltonians at all. In contrast to it, 
in the parafermionic approach on the one hand the type A 3-fold superalgebra is always 
realized in the subsector with the fermion number and 3, and on the other hand the vari- 
ous patterns in the existence of intermediate Hamiltonians are characterized by considering 
the other subsector with the fermion number 1 and/or 2. In addition, generalized M- fold 
superalgebra like (7.8) could provide an alternative for defining paraSUSY of order Af. In 
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Case Class (1, 1) 



Class (0,1), (1,0) 



I a/2 

II 2z 

ir 

III z 2 /2 

IV ±2a(z 2 - 1) 



2z 3 



±2a(z 2 + 1) 

IV 



2az 2 (z + 1) 

2z{l - z)(l - mz) 

z{z 2 +2(l-2m)z + l)/2 



V 



TABLE III: A real classification scheme of type A 3-fold SUSY models with one set of intermediate 
Hamiltonians. In the above, a 6 R is an arbitrary constant and < m < 1. 

the present A/" = 3 case, the conventional defining algebra (7.1b) only produces in essence 
the additional constraint (7.6) whose physical relevance is unclear. Hence, there seems, at 
least until now, no physical evidence to claim that which definition is appropriate. 

In principle, one can continue to study the Af > 3 cases, but as already mentioned the 
number of different patterns in the existence of intermediate Hamiltonians drastically in- 
creases as Af increases. However, the number of different sets of intermediate Hamiltonians 
would be limited to at most 4 regardless of the number of patterns due to the following 
reasons. First, the transformation formula (2.10) which is valid for all Af > 3 could produce 
algebraic equations of at most fourth-degree and thus at most four additional sets of interme- 
diate Hamiltonians would be admissible. But type A AAfold SUSY systems for Af > 3 have 
at most three independent free parameters bi (i = 0, 1,2) and thus at most three different 
solutions among the four would be compatible simultaneously. Hence, at most three addi- 
tional sets would be available, which means that the maximum number of different sets is 
four in total. By a similar argument we conclude that it would be at most 3 if one constraint 
on the parameters 6, are inevitable for the existence of one set of intermediate Hamiltonians, 
as in the present Af = 3 case (4.1)-(4.3), since in this case there are essentially at most two 
independent free parameters and thus at most two algebraic solutions would be compatible 
simultaneously. The latter fact is indeed the reason why in the Af = 3 case there are at most 
three different sets and thus are no classes such as Class (to, m + n) with to + n > 3, Class 
(to + I, to + n) with to + I + n > 3, Class (0, n) with n > 3, and Class (I; n) with I + n > 3. 
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